Observation of Long-Lived Muonic Hydrogen in the 2S State 
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The kinetic energy distribution of ground state muonic hydrogen atoms /ip(lS') is determined 
from time-of-flight spectra measured at 4, 16, and 64 hPa Eb room-temperature gas. A 0.9 keV- 
component is discovered and attributed to radiationless deexcitation of long-lived fip(2S) atoms in 
collisions with Hb molecules. The analysis reveals a relative population of about 1 %, and a pressure- 
dependent lifetime (e.g. 30.4 Ig 1 ? 4 ns at 64 hPa) of the long-lived fj,p(2S) population, equivalent to 
a 25'-quench rate in (ip(2S)+H.2 collisions of 4.4 j^'g x 10 11 s _1 at liquid hydrogen density. 

PACS numbers: 36.10.Dr, 34.20.Gj 



A measurement of the Lamb shift in muonic hydrogen, 
i.e. the energy difference of 0.2 eV between the 2P- and 
2S'-states of /i - p atoms 0, 0, is m progress at the 
Paul Scherrer Institute (PSI), Switzerland 4]. Vacuum 
polarization shifts the 25'-levels by 0.2 eV below the 2P- 
levels; fine and hyperfine splittings of the n = 2 levels are 
much smaller. The finite size effect is 2 % of the Lamb 
shift. 

The /ip Lamb shift experiment is expected to give a 
precise value for the root-mean-square charge radius of 
the proton [j| . Together with recent advances in H-atom 
spectroscopy |6J, [7| , this leads to a better determination 
of the Rydberg constant, and to a test of bound-state 
quantum electrodynamics on a new level of precision Q . 

The most important prerequisite for such an exper- 
iment, the availability of sufficiently long-lived /ip(25) 
atoms, has so far not been experimentally established. 
When muons are stopped in H 2 gas, /ip atoms are formed 
at high n-levels and then deexcite predominantly to the 
ground state ("muonic cascade"). A fraction £25 of a 
few percent reaches the metastable 2S state whose life- 
time is, in absence of collisions, essentially given by the 
muon lifetime of 2.2 /is. 

In a gas, there is collisional 2S'-quenching, with very 
different rates depending on the /xp kinetic energy be- 
ing above or below the 2S-2P threshold («0.3 eV in the 
lab frame) ||. Most uv(2S) atoms are formed at ener- 
gies above this threshold |l0j , where collisional 25 — -> 2P 
Stark transitions (followed by 2P — *1S radiative deex- 
citation) lead to rather fast 2 S- depletion. This is the 
"short-lived" 2S-cornponeiit with a predicted lifetime flU 
r 2S° rt ~ I00ns/pH 2 [ n Pa]- There is, however, a competi- 
tion between such Stark transitions and deceleration 0. 
A fraction of the fip(2S) atoms should therefore survive 
the process of slowing down below 0.3 eV, where transi- 
tions to the 2P state are energetically forbidden. These 
Hp(2S) atoms form the "long-lived" 2S'-component, with 
a lifetime T2™ s and a population e^™ 8 (P er A*P atom). 



The 25-population £25 is well determined from the 
measured ^p X-ray yields 0,EJ0|- el^can be derived 
from the measured IS* kinetic energy distributions by 
using calculated elastic and inelastic /ip(25) cross sec- 
tions 0. At 16 hPa, for example, e 2 s = (4.40 ± 0.17) % 
and £ 2 °s S = (F16 ± 0.12) %. Calculations of the radiative 
quenching process during collisions |16L Il7l Il8| predicted 
values of t^™ 8 ~ 600 ns or larger for an H2 pressure of 
64hPa, but this could not be confirmed e xper imentally 
in searches for delayed n = 2^1 X-rays [12J, |l_9 , 20] . 

We report here on the first measurement of a non- 
radiative deexcitation mechanism for long-lived fip(2S) 
atoms at gas pressures ph 2 = 4, 16, and 64 hPa at room 
temperature |l5| . It was found in the analysis of time-of- 
flight (TOF) spectra which were measured to determine 
the /ip(15) kinetic energy distributions at low densities. 
A pronounced fast component showed up, corresponding 
to a surprisingly high /ip(lS') kinetic energy of 0.9 keV. 

The basic idea of the experiment is to form /xp(lS') 
atoms near the axis of a cylindrical gas target and to de- 
duce their kinetic energy from the measured TOF they 
need to reach the gold-coated cylinder surface. The muon 
transfer reaction /ip +Au — > /iAu*+p to highly excited 
/iAu* takes place, and several /iAu X-rays with energies 
up to 9 MeV are emitted and detected by a CsI(Tl) scin- 
tillator. The cylindrical geometry gives a better kinetic 
energy resolution than earlier planar geometries [2l| . 

Figure ^ shows the experimental setup. Negative 
muons with a momentum of p = 10 MeV/c from the high- 
intensity 7rE5 channel at PSI entered aim long, 20 cm 
inner diameter solenoid operated at 5 T and evacuated 
to 10 -5 hPa. Its fringe field reduced the radial beam 
size to a few mm in the region of full field. About 
5000 /i~/s reached a 16 /jm thick Mylar entrance mod- 
erator. The muons were detected by a parallel-plate 
avalanche counter (PPAC) with a total areal density of 
300 /ig/cm 2 . The PPAC is a double gas chamber filled 
with 16 hPa isobutane, with a detection efficiency up to 
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FIG. 1: Experimental setup inside the 5 T solenoid (not 
to scale, dimensions given in text). PPAC: parallel plate 
avalanche counter; PS: plastic scintillator; MCP: microchan- 
nel plate; PMT: photomultiplier tube. U(z) is the electric 
potential of the 23 stack rings (7 with C-foils). 



80%, depending on the muon energy. At the exit of the 
PPAC, the muons had energies between and 200 keV. 
They then crossed a stack of seven thin carbon foils (each 
with an areal density of 5 fig/ 'cm 2 and 2 cm diameter) 
where a few electrons were released. A high voltage ap- 
plied to the stack (AU = 1800 V between successive foils) 
accelerated these electrons, and frictional cooling [22| en- 
hanced the muon stop rate. After traversing the gas tar- 
get, the electrons were detected in a microchannel plate 
(MCP). An incoming muon was identified by a delayed 
coincidence between the signals from PPAC and MCP. 
The muon detection efficiency of the stack-MCP assem- 
bly was ~ 40 % for a target gas pressure of 4 hPa and 
decreased to ~ 15 % at 64 hPa. The radial beam distribu- 
tion was measured to be approximately Gaussian-shaped 
with spreads a x ~ <r y ~ 1.8 mm in the 5 T field. The 
cyclotron radii of the muons were ~ 1 mm at 5 T. 

The gas target vessel was a 26 cm long Al cylinder. Its 
inner cylindrical surface was coated with a 500 nm thick 
evaporated gold layer. The side walls were covered with 
polypropylene. Measurements were performed with two 
target inner radii, R — 29 mm and 10 mm. A target with 
polypropylene instead of gold was used to determine the 
background. Muon transfer to carbon does not produce 
X-rays within the accepted /^Au X-ray energy window. 

The upstream and downstream target windows were 
thin Formvar foils (1 cm diameter) of 10 /xg/cm 2 at pn 2 = 
4 and 16 hPa, and 40/Ltg/cm 2 at 64 hPa. H2 gas with 
< 3 ppm impurities was continuously flushed through the 
target (exchange rate 2/min). Out gassing from walls and 
pipes resulted in a total impurity concentration of order 
10 -4 . The measured /ip formation rate was ~20 s _1 . 

A 20 cm long cylindrical CsI(Tl) scintillator of 9.2 cm 
inner and 19 cm outer diameter was placed around the 
target. It had a high probability to detect at least one of 
the /iAu X-rays. The CsI(Tl) crystal was surrounded on 
all but the outer side by an 8 mm thick plastic scintillator 
acting as anti-coincidence detector against charged par- 
ticle background. The fast light component of the plastic 
scintillator was distinguished from the "slow" Csl light 
component (decay time ~ 1.5 /js) by pulse shape analysis. 



For each incoming muon, the time difference between 
the MCP signal and a Csl signal from a fj,Au X-ray was 
recorded. The true /zp(lS) TOF, however, starts at the 
formation of a np(lS). Time calibration measurements 
were done with a few per cent of xenon admixed to the H2 
target gas, because at such a low Xe concentration /iXe 
and /ip formation times are nearly identical. The /xXe X- 
rays were detected in the CsI(Tl) crystal, too. Several Xe 
concentrations were used and the fip(lS) formation time 
distribution was obtained by extrapolating to zero Xe 
concentration. A careful analysis of the /xXe time spectra 
separated the time distribution of direct /zXe-formation 
from the (delayed) transfer reaction /xp +Xe — > /xXe*+p. 
Small residual deviations not taken into account by the 
Xe calibration (like the /ip cascade time and discrimina- 
tor threshold effects) were also considered. 

Figure [3 shows the fiAu time spectra measured at 64, 
16, and 4 hPa for both target radii, and the ^tp for- 
mation time distributions deduced from the /iXe data. 
The background was studied with the polypropylene tar- 
get. It is dominated by bremsstrahlung from muon-decay 
electrons, with an additional component (depending on 
pressure) at times preceding /xp formation, correlated to 
muons stopping in high-Z materials near the target. All 
/iAu time spectra exhibit a wide distribution whose maxi- 
mum corresponds to /xp kinetic energies Skin of 5... 20 eV. 
Such energies result from Coulomb deexcitations [2^ 
/xpn+Ha — » MPn'<n + ^kin +P + ••• during the cascade, a 
well known effect for /ip and 7rp atoms |9j. 

The data measured in the R = 29 mm target, which 
is most sensitive to high kinetic energies, show an un- 
expected peculiarity: The 64 hPa spectrum features an 
additional pronounced peak at t w 100 ns. This peak is 
also visible at 16 hPa, but almost disappears at 4 hPa. 
The corresponding kinetic energy is about 0.9 keV, in- 
dicated in Fig. Elas dotted line at 0.074 fj,s which is the 
calculated TOF for 29 mm, the shortest distance from 
the beam center to the Au wall. 

Such a high-energy /j,p(lS')-component has never been 
reported before. It can only originate from radiation- 
less deexcitation of /Ltp(n = 2) atoms in collisions with H2 
molecules. The n = 2 — » 1 transition energy of 1.9 keV is 
shared according to momentum conservation between the 
formed /ip(lS) and a proton of the H 2 molecule: 

/ip(n = 2) + H 2 -> W is (0.9 keV) + p(l keV) + . . . . 

A detailed analysis of the measured time spectra was 
performed using a set of Monte Carlo (MC) generated 
TOF spectra extending over the whole expected energy 
range (~ 10~ 2 - 10 3 eV). In the MC, the target and detec- 
tor geometries, axial and radial muon stop distributions, 
X-ray detection efficiencies, and scattering of fip atoms 
on H2 molecules and on the gold surface of the target 
cylinder were taken into account (l5| . Cross sections for 
/xp(15 , )+H2 scattering calculated by Adamczak |24| for 
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FIG. 2: (color online). Early part of the /uAu time spec- 
tra measured at 64, 16, and 4 hPa (top to bottom) with 
the R = 29 mm (left) and i?=10mm (right) targets. The 
corresponding pip kinetic energies are also indicated. Com- 
ponent (a) at 0.9 keV originates from non-radiative pp(2S)- 
quenching and was convoluted with the fitted 2<S'-lifetime. 
The thin solid line is the total fit function, composed of the 
measured background (thick black line), the 2S component 
(a), the high energy Coulomb deexcitation lines with main 
components (b) n — 3 — » 2 (166 eV) and (c) n = 4 — > 3 (58 eV), 
and a continuous distribution at energies below 60 eV (d). 
The dashed curve (s) is the fj,p formation time distribution 
obtained in the Xe calibration measurements. 



gent consistency test because the two time spectra are 
very different for the same kinetic energy distribution. 

An adequate fit of the data could not be obtained un- 
less it was assumed that the 0.9 keV component was 
emitted with a non-zero lifetime. Such an effect can 
only be caused by long-lived /xp(25) atoms. (Other colli- 
sional processes, like muon transfer to medium- Z atoms 
from gas impurities, would produce X-rays with ener- 
gies below those accepted for fxAu.) For the data anal- 
ysis a two-step process was adopted where the fip(2S) 
atoms are first slowed down to low energies without sig- 
nificant quenching, and then collisionally quenched at a 
constant rate Ajlg 6110 . The 0.9 keV MC time spectrum 
was therefore convoluted with an exponential, with life- 



time r. 



long 
25' 



(A' 



quench 
2S 



+ Xfj.) 1 and amplitude A 2 s 



£ 2°s S A2s° nch ; delayed by a time offset At 2 $ for slowing 
down. (A M =4.55 x 10 5 s _1 is the muon decay rate.) 

First the data sets at 64 and 16 hPa were analyzed 
independently of each other. (At 4 hPa the 0.9 keV 
component is too small to be treated in an indepen- 
dent fit.) Free fit parameters were the 25-quench time 
T 2S° nCh = l/A^ 011011 , t ue amplitude A 2 $, the time offset 
A^25, and the amplitudes of the Coulomb deexcitation 
components. The results are given in Tab. |H For At 2S 
a value of 8 ns was obtained at 16 hPa, and roughly 
4 times less at 64 hPa, i.e. much smaller than 7"2<j ench . 

Next, T2s Cnch was set to be proportional to l/pn 2 , be- 
cause non-linear effects (like three-body collisions) can be 
excluded at our conditions. A combined fit of the data 
at all three pressures gave nearly the same values as the 
independent fit (Tab.|IJ. The corresponding quench rate 
at room temperature is 



A 



quench 
2S 



5.1 



+2.4 
-2.1 



x 10 5 



Ph 2 [hPa] 



or, normalized to liquid hydrogen atom density 
(LHD = 4.25 x 10 22 atoms/cm 3 ), 



A~ h (LHD) = AAtf.l x 1Ql1 s_1 - 



(1) 



laboratory kinetic energies below 9 eV («s H2 breakup en- 
ergy) have been used. Between 9 eV and 1000 eV we ap- 
plied those calculated by Cohen |2f| for collisions with 
H-atoms. Coulomb deexcitation lines from low n-levels 
were included in the fit function; the most important 
transitions are indicated in Fig. [21 A free fit gave only 
small amplitudes for Coulomb deexcitation components 
with An > 1, in accordance with calculations |23|. Fi- 
nally, each MC TOF spectrum was convoluted with the 
measured fip formation time distribution. The fit func- 
tion is a superposition of these spectra, with their in- 
dividual weights treated as free parameters. For each 
pressure, the two time spectra measured for R = 10 mm 
and 29 mm were fitted simultaneously, providing a strin- 



2S'-quench time T^ cnch : 



-1 / \ quench 
: V A 2S 



and 



TABLE I: 

population £2°s g (per /ip atom) of the long-lived fip(2S)- 
component, obtained for the "independent" and "combined" 
fits of the time spectra. For the combined fit we give the 
statistical and the systematic errors, respectively. Only the 
statistical errors are given for the independent fits. 
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combined fit 
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+0.69 +0.29 
-0.25 -0.23 
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The errors result from the quadratic sum of statisti- 
cal and systematic errors. The combined fit gave a re- 
duced x 2 of 1.030 (x 2 = 529.5 for 514 degrees of freedom). 
X 2 = 559.2 resulted for T^ cnch (64 hPa) fixed to 9 ns, i.e. 
such a short 2S'-lifetime is statistically excluded at the 
5er confidence level. 

Various sources of systematic errors were investigated: 
/zp formation time, scattering on H2 and gold, radial size 
of the muon stop distribution, and background shape. 
Uncertainties in the formation time gave a contribution of 
ns to r25° nch (64 hPa), H 2 -scattering one of ns, 
whereas all other errors were considerably smaller. Gas 
impurities modified T2g Cnch (64 hPa) by < 0.1 ns. The re- 
sulting total systematic errors are given in Tab. [I] 

The large value found for A^ 011 ^ can be attributed to 
the formation of an excited muonic molecular ion in a 
resonantly enhanced reaction with H2 jl^, [2]j 

/i P (2S)+H 2 ^{[(pp M )+]*pee}* , 

similar to the well-known formation of ground state 
muonic molecules in muon catalyzed fusion |28j| . There 
is a series of resonant states just below the fJ-p(n = 2) + p 
dissociation limit |29j; rovibrational excitations of the 
electronic molecule absorb the small excess energy from 
the formation of the (pp^t) + "pseudonucleus" . A subse- 
quent Auger decay to states with > 16 eV binding energy 
stabilizes the [(pp//) + ]* against back decay. Mixing of the 
excited [(pp/i) + ]* with the repulsive (pp/i) + ground state 
finally leads to /ip(lS) atoms with 0.9 keV kinetic energy. 

A theoretical estimate of the molecular formation 
rate [3(| gave ~5 x 10 10 s" 1 at LHD, ten times smaller 
than our result (Eq. Q). The calculation neglected the 
correlated motion of the three interacting atoms which 
is important in low-energy fip(2S) collisions because 
of its long-range polarization potential. The triatomic 
molecule has a much denser vibrational-state spectrum 
with more ways to exchange excitation energy, lea ding 
to longer effective lifetimes of the compound system |3jJ 
and correspondingly higher molecular formation rates. 

The populations e 2 °s S resulting at 16 and 4 hPa agree 
with the values of (1.16 ± 0.12)% and (1.03 ± 0.09)%, 
respectively, evaluated for the /j,p(2S) that slow down 
to below 0.3 eV 0|. This means that there is little 
room for radiative 2S — > IS deexcitation processes and 
explains why no delayed 1.9 keV X-rays were found in 
previous experiments |l2Hl9Ll2fij |. The dominance of non- 
radiative 2«S'-decay resulting from our measurements is in 
accord with recent calculations of the radiative and non- 
radiative decay rates of [(pp/x) + ]* molecules |29ll32]| . 

In summary, we have measured that ~ 1 % of all 
muons stopped in H2 gas between 4 and 64 hPa form 
long-lived fip(2S) atoms that decay non-radiatively into 
fip(lS) atoms of 0.9 keV kinetic energy. For the cur- 
rent /Ltp Lamb shift experiment, a pressure of ~ 1 hPa 
is needed where both 2S'-population (~ 1 %) and lifetime 



( T 2s" S = 1-04 -0 21 M s & t 1 hPa) meet the experimental re- 
quirements 0. First measurements at a new low-energy 
muon beam at PSI showed that enough /zp atoms can be 
formed in a small gas volume at such pressures. 
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